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Abstract — The spectrum of the Laplacian matrix of a network 
plays a key role in a wide range of dynamical problems associated 
with the network, from transient stability analysis of power 
networks to distributed control of formations. Using methods 
from algebraic graph theory and convex optimization, we study 
the relationship between structural features of a network and 
spectral properties of its Laplacian matrix. We illustrate our 
results by studying the influence of structural properties on 
the Laplacian eigenvalues of the American (western states), 
French and Spanish high-voltage transmission networks. Our 
study suggests that for such networks the Laplacian spectral radii 
and spectral moments are strongly constrained by a particular 
set of local structural features, namely the degree sequence and 
the so-called joint-degree distribution. On the other hand, other 
structural properties that may seem important, such as the 
distribution of cycles in the network, appear to have a very weak 
influence on the Laplacian spectrum of electrical transmission 
networks. We also show that local structural features are not 
enough to characterize the Laplacian spectral gap. Therefore, 
since the spectral gap is fundamental in the analysis of many 
dynamical processes on networks, random models in which only 
local structural features are prescribed are typically insufficient 
to generate synthetic topologies in which these dynamical pro- 
cesses can be studied. 



I. Introduction 

The eigenvalue spectrum of the Laplacian matrix of a 
network provides valuable information regarding the behavior 
of many dynamical processes taking place on the network. 
In their seminal paper ([Tl, Pecora and Carrol related the 
problem of synchronization in a network of coupled oscillators 
to the largest and second-smallest Laplacian eigenvalues of 
the network (usually denoted by Laplacian spectral radius and 
spectral gap, respectively). More recently, Dorfler and BuUo 
derived conditions for transient stability in power networks in 
terms of the spectral gap of the Laplacian matrix Q. Apart 
from their applicability to the problems of synchronization and 
transient stability analysis, the Laplacian eigenvalues are also 
relevant in the analysis of many distributed estimation and 
control problems (see fSl,!?! and references therein). 

Understanding the relationship between the structure of a 
complex network and the behavior of dynamical processes 
taking place in it is a central question in the research field of 
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network science [5|. Since the behavior of many networked 
dynamical processes is closely related with the Laplacian 
eigenvalues ||6l, it is of interest to study the relationship 
between structural features of the network and its Laplacian 
eigenvalues. In this paper we study this relationship, focusing 
on the role played by structural features that can be extracted 
from localized samples of the network structure. Our objective 
is then to efficiently aggregate these local samples of the 
network structure to infer global properties of the Laplacian 
spectrum. Our analysis reveals that there are certain spectral 
properties, such as the spectral radius and the so-called spectral 
moments, that can be efficiently estimated from local structural 
features of the network. 

The most common approach to study the relationship be- 
tween structural and spectral properties of a network is by 
means of random network models. This approach has been 
widely used to study the effect of the degree distribution |7|, 
correlations O, or clustering ||9l, in a broad range of real- 
world networks. Although very common in the literature, this 
approach presents a serious drawback: when we prescribe the 
structural feature under study in the random network model, 
there are other structural features that are being indirectly 
induced in the network and may be relevant in the network's 
behavior Since those induced features are not directly con- 
trolled, it is difficult (if not impossible) to isolate the role 
of the structural feature under study using random network 
models. 

In this paper, we present a novel approach to study the 
influence of certain structural features on the network Lapla- 
cian eigenvalues without using random models. Our approach 
builds on algebraic graph theory and convex optimization 
to study relevant properties of the Laplacian spectrum. We 
apply our theoretical results in the study of real electrical 
transmission networks. In particular, we analyze high- voltage 
transmission networks in the USA (western states), France 
and Spain. Our numerical results show that, in the case of 
electrical transmission networks, the Laplacian spectral radii 
and spectral moments are strongly constrained by a particular 
set of local structural features: the degree distribution and 
the so-called joint-degree distribution. On the other hand, 
other structural properties that may seem important, such as 
the distribution of cycles in the network, appear to have a 
very weak influence on the Laplacian spectrum of electrical 
transmission networks. In our numerical experiments, we also 
verify that local structural features are not enough to charac- 
terize the Laplacian spectral gap, since this spectral property 
measures, to a certain extent, how well a network is globally 
connected. Therefore, since the spectral gap is fundamental 



in the analysis of many dynamical processes on networks, 
random models in which only local structural features are 
prescribed, such as those in Q-lEI, are insufficient to generate 
synthetic topologies in which these dynamical processes can 
be studied. 

The rest of this paper is organized as follows. In the 
next section, we define graph-theoretical terminology needed 
in our derivations. In Section |lllj we review some previous 
results relating the eigenvalues of the Laplacian matrix of 
a network to its structural features, and introduce a novel 
methodology based on algebraic graph theory to derive closed- 
form expressions for the so-called Laplacian spectral moments. 
In Section |IV| we use convex optimization to derive tight 
bounds on important properties of the Laplacian spectrum 
from the spectral moments. We apply our theoretical results 
to the analysis of the Laplacian spectra of the American 
(western states), French and Spanish high-voltage transmission 
networks in Section [Vl 

II. Preliminaries 

Let Q = {V,£) be an undirected graph, where V = 
{vi, . . . ,Vn} denotes a set of n nodes and £ C V x V denotes 
a set of e undirected edges. If {vi,Vj} £ £, we call nodes 
Vi and Vj adjacent (or first-neighbors), which we denote by 
Vi ^ Vj. We define the set of first-neighbors of a node Vi as 
Afi ~ {w E V : {vi,w} E £}. The degree di of a vertex 
Vi is the number of nodes adjacent to it, i.e., di = \JVi\. We 
consider three types of undirected graphs: 

(/) A graph is called simple if its edges are unweighted 
and it has no self-loopaH 

(;/) A graph is loopy if it has self-loops. 

(Hi) A graph is weighted if there is a real number asso- 
ciated with every edge in the graph. 

More formally, a weighted graph H can be defined as the 
triad H = {V,£, W), where V and £ are the sets of nodes and 
edges in H, and W = {wij e M, for all {vi,Vj} e £} is the 
set of (possibly negative) weights. 

Graphs can be algebraically represented via matrices. The 
adjacency matrix of a simple graph Q, denoted by Ag — [uij], 
is an n X n symmetric matrix defined entry-wise as a,j- — 1 
if nodes Vi and Vj are adjacent, and a^ = otherwise. 
In the case of weighted graphs (and possibly non-simple), 
the weighted adjacency matrix is defined by Wg ~ [wij], 
where Wij = if Vi is not adjacent to Vj. We define the 
degree matrix of a simple graph Q as the diagonal matrix 
Dg — diagidi). We define the Laplacian matrix Lg (also 
known as combinatorial Laplacian, or Kirchhoff matrix) of 
a simple graph as Lg ~ Dg — Ag. For simple graphs, 
Lg is a symmetric, positive semidefinite matrix, which we 
denote by Lg ^ ifTOl . Thus, Lg has a full set of n real 
and orthogonal eigenvectors with real nonnegative eigenvalues 
= Ai < A2 < ... < An. Furthermore, the trivial eigenvalue 
Ai = of Lg always admits a corresponding eigenvector 
Vl = (1,1,...,!) . The algebraic multiplicity of the trivial 
eigenvalue is equal to the number of connected components 

'a self-loop is an edge of the type {vi, Vi}. 



in Q, which we assume to be equal to one in the rest of the 
paper The smallest and largest nontrivial eigenvalues of Lg, 
A2 and A„, are called the spectral gap and spectral radius 
of Lg, respectively, and they play an important role in this 
paper Given a nx n real and symmetric matrix B with (real) 
eigenvalues ci, ..., tT„, we define the fc-th spectral moment of 
B as. 



1 

n ^ — ^ 



As we shall show in Section III-B there is an interesting 
connection between the spectral moments of the Laplacian 
matrix, mk [Lg), and structural features of the graph. 

We now define a collection of structural properties that are 
important in our derivations. The degree sequence of a simple 
graph Q is the ordered list of its degrees, (di, ..., (i„). The p-th 
power-sum of the degree sequence is defined as 
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A walk of length k from vi^ to Vi^^^ is an ordered sequence 



of nodes [vi^ , vi 



..,v 



tk+i, 



such that Vi 



Vtj+i for j = 



1, 2, ..., fc. One says that the walk touches each of the nodes 
that comprises it. If w^^ — Wi^^^, then the walk is closed. 
A closed walk with no repeated nodes (with the exception 
of the first and last nodes) is called a cycle. Given a walk 
p = (vi-^ , Vi^, ..., Wifc_|_i) in a weighted graph H, we define the 
weight of the walk as, uj (p) = Wi^i^Wi^i^...Wi^i^^^. 

Finally, we define the concept of local neighborhood around 
a node. Let 5 {vi, Vj) denote the distance between two nodes 
Vi and Vj (i.e., the minimum length of a walk from Vi to Vj). 
We say that Vi and Vj are r-th neighbors if 5{vi,Vj) — r, 
and define the r-th order neighborhood around Vi as the set 
of nodes Af^ ^ {w E V : 6 {vi,w) < r}. By convention, we 
assume that S{v,v) ~ 0, hence, the first-order neighborhood 
of Vi satisfies A/"/ = AfiU Vi. The nodes in A/'j induce an 
(unlabeled) local subgraph Qy"' ~ (A/'j , £j ) C Q, where 
^i — {{^1 w} G £ s.t. v,w E Mi }, i.e., the set of edges 
connecting nodes in A/^ . It is worth remarking that these 
local subgraphs are unlabeled; thus, it is impossible, in general, 
to reconstruct the complete network Q from the set of local 
subgraphs {Gi , * = 1, ...jJt-}. These subgraphs are useful to 
define the concept of local structural measurement. We say 
that a structural measurement is local with a certain radius r 
if it can be computed from the set of local subgraphs {Ql , 
i ~ \, ..., n}. For example, the degree sequence of a graph is 
a local structural measurement (with radius 1), since we can 
compute the degree of each node Vi from the local subgraph 
A^^ ' . Similarly, one can compute the number of triangles 
touching Vi from M^ ; hence, the total number of triangles in 
t/ is a local measurement with radius 1. In contrast, the eigen- 
value spectrum of the Laplacian matrix is not a local property, 
since we cannot compute the eigenvalues unless we know the 
complete graph structure. One of the main contributions of 
this paper is to propose a novel methodology to extract global 
information regarding the Laplacian eigenvalue spectrum from 
local structural measurements of the network. 



III. Spectral Analysis of the Laplacian Matrix 

The spectrum of the Laplacian matrix is relevant in a wide 
range of networked dynamical processes, such as synchro- 
nization of coupled oscillators |1|, flocking and formation 
control of multi-agent systems E),!!], or distributed consensus 
algorithms fTTl. In the context of power networks, Laplacian 
eigenvalues contain useful information for the analysis of 
synchronization and transient stability llj. Among the Lapla- 
cian eigenvalues, the spectral radius and the spectral gap are 
specially important in dynamical applications and we pay 
particular attention to them. 

In this paper we introduce a novel approach to study the 
effect of structural properties of a network on the spectrum of 
its Laplacian matrix, by using the moments of the eigenvalue 
distribution. Before we introduce our approach, we first review 
some of the most relevant known bounds on the Laplacian 
eigenvalues. 



As we shall verify in Section [Vj the above upper bound on 
the spectral radius is quite tight in the case of electrical trans- 
mission networks, while lower bounds found in the literature 



are very loose. In Section IV we shall derive lower bounds 
that outperform those found in the literature in the case of 
electrical transmission networks. 



B. Moment-Based Analysis of the Laplacian Spectrum 

In this section, we use algebraic graph theory to relate 
spectral properties of a network to a rich variety of local 
structural measurements. First, we need to introduce some 
preliminary results and concepts. A well-known result in 
algebraic graph theory relates the diagonal entries of the k- 
th power of the adjacency matrix, [A^] .., to the number of 
closed walks of length k m Q that start and finish at node 
Vi [10|. We can generalize this result to weighted graphs as 
follows: 



A. Bounds on Laplacian Eigenvalues 

The literature contains a large collection of upper and lower 
bounds on the spectral gap A2 and the spectral radius A„ 
in terms of structural features of a network, such as degree 
distributions. Let us consider a simple graph Q with n nodes, 
e edges, and degree sequence {di)^^^. Let m^ be the average 
degree of the neighbors of w^, i.e., m^ = J" '^v-eM- '^i (note 
that rrii is a local structural measurement with radius 2). 
For connected networks, we have the following upper and 
lower bounds on the Laplacian spectral gap in terms of local 
structural measurements (see [12,1 for an extensive review): 
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(n — 1) m 
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A2 < Ke < min {di} , from lfT4l . 
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where ra = n'^v-ev^^'- ^'^'^ '^e ^^ '■^^ edge connectivity^ 
of the graph. Although many other bounds can be found in 
the literature, most of them involve properties that cannot 
be computed from local structural features of the network, 
such as the diameter of the network or the average distance 
between every pair of nodes (12]. Since the spectral gap 
strongly depends on the global connectivity of the network, 
bounds on the spectral gap written in terms of local structural 
features usually perform poorly in practice (as we shall verify 
in Section \V\. 

Similarly, we find in the Uterature a large variety of upper 
and lower bounds on the Laplacian spectral radius A„ (see, 
for example, |15| for a collection of bounds). Many of these 
bounds can be written as the maximum or minimum of certain 
functions defined over the set of nodes or edges. For example. 



di + dj + y {di — dj) + Arriimj 
A„ < max , |[T6l, 

{vi,Vj}££ A 



^n ^ V 2 min (dj + diT 



>l/2 
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Proposition L- Let 7i={V,£,W) be a weighted graph 
with weighted adjacency matrix W-h = [wij]. Then 



where P^;*^ is the set of closed walks of length k from Vi to 
itself in the complete loopy graplrl and u (p) is the weight of 
walk p in H. 

Proof: By recursively applying the multiplication rule for 
matrices, we have the following expansion 



M„- 



n n 

EE- 

i=l 12 = 1 



■E^ 



(2) 



Using the graph-theoretic nomenclature introduced in Sec- 
tion III] we have that Wi^i^Wi^ i^.-.Wi^j = w{p), for p — 
{vi,Vi^,Vi^, ...,Vi^,Vi). Hence, the summations in S can be 
written as [W^]^. = J2l<i,^2,....^,<n'^ (p)- Finally, the set 
of closed walks p = {vi,Vi2,Vi^, ...,Vi^,Vi) with indices 
1 < i, 12, •■•, ik < nis equal to the set of closed walks of length 
k from Vi to itself in the complete loopy graph J„ (which we 
have denoted by P^'^ in the statement of the Proposition). ■ 

The above Proposition allows us to write the spectral 
moments of the weighted adjacency matrix of a weighted 
graph H in terms of the weights of closed walks, as follows: 

Lemma 3.1: Let Ti. = {V,£,W) be a weighted graph with 
weighted adjacency matrix W-h = [wij]. Then 



^k{Wn)^-y2 E ^b): 



V,€Vp^p(^)^ 



where P^: ^ is the set of closed walks of length k from Vi to 
itself in the complete loopy graph J„. 



-The edge connectivity of a connected graph is the minimum number of ^The complete loopy graph Jn is the graph with node-set V = {vi, ..,,Vn} 

edges whose removal renders a disconnected graph. and edge-set £ = {{vi,Vj} for all Vi,Vj G V}. 
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Fig. 1. A simple graph 54 (left) and its corresponding Laplacian weighted 
graph C{Qi) (right). 



Proof: Let us denote by /ii,...,/i„ the set of (real) 
eigenvalues of the (symmetric) weighted adjacency matrix 
W-u- We have that the moments can be written as 



TOfc 



1 " 1 

(M/«) = -Em'- -Trace (ly^) 



i=l 



since W-u is a symmetric (and diagonalizable) matrix. We then 
apply Proposition [T] to rewrite the moments as follows. 



rrik 



(^«) = ^E[<]. = ^E E -w- 
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In Subsections III-Bl and III-B2 we shall apply this result 
to compute the spectral moments of the Laplacian matrix 
in terms of local structural measurements. First, we need to 
introduce a weighted graph that is useful in our derivations: 

Definition 3.1: Given a simple graph Q ={y,£), we de- 
fine the Laplacian graph of Q as the weighted graph 
C [g) = (V, 5 U SnX), where 5„ = {{w, v} for all v G V} 
(the set of all self-loops), and F = [7^^] is a set of weights 
defined as: 



Irj 



— 1, for {vi,Vj} E £ 
di, for i=j 
0, otherwise. 



Remark 3.1: Note that the weighted adjacency matrix of the 
Laplacian graph C{Q) is equal to the Laplacian matrix of the 
simple graph Q. Hence, we can apply Lemma 3.1 to express 



the spectral moments of the Laplacian matrix Lg in terms of 
weighted walks in the Laplacian graph C {Q). 

Example 3.1: In Fig. [TJa) we have a simple graph with 4 
nodes, Q/^. In Fig. flTb) we have its corresponding Laplacian 
weighted graph C [Ga)- Observe how the edges in Q^ are also 
edges in C {Q4) with an associated weight of —1. Also, each 
node in C {Q4) presents a self-loop with a weight equal to 
the degree of the corresponding node in CJ4. One can easily 
check that the Laplacian matrix of Q/^ is equal to the weighted 
adjacency matrix of C {Gi)- 

Before we apply Lemma [TT] to study the Laplacian spectral 
moments, we must introduce the concept of subgraph covered 
by a walk. 

Definition 3.2: Consider a walk p ~ (vi^ ^Vi^ Wi,, , , ) of 



"■^k+i) 



length fc in a (possibly loopy) graph. We define the subgraph 



covered by p as the simple graph C [p] = (Vc [p) ,£c (p)), 
with node-set Vc (p) 



Ur=i '"ir^ ^^'^ edge-set 



Sc (P) = 



u 



b 



Hr + 



J , for 1 < r < fc + 1. 



Based on the above, we define triangles, quadrangles and 
pentagons as the subgraphs covered by cycles of length three, 
four, and five, respectively. Notice that self-loops are excluded 



from £c {p) in Definition 3.2 For example, consider a walk 
P = (vi,vi,V2,v-j„V'i,Vi,v-i,vi) in a graph with self-loops. 
Then, C (p) has node-set Vc (p) — {vx^v-i^v-^ and edge-set 
^c (p) = {{v\,vi\ , {w2, W3} , {u3, wi}}. In other words, C (p) 
is a simple triangle. 

to 



In the following subsections, we build on Lemma 3.1 



derive closed-form expressions for a sequence of Laplacian 
spectral moments in terms of relevant structural features of the 



network. In Subsection III-B 1 we provide expressions for the 
first three spectral moments in terms of the degree sequence 
and the number of triangles in Q. In Subsection III-B2 we 



compute higher-order moments in terms of longer cycles and 
more elaborate structural measurements. 

1) Low-Order Laplacian Spectral Moments: In this subsec- 
tion we derive expressions for the first three spectral moments 
of the Laplacian matrix of a simple graph using graph- 
theoretical concepts. 

Theorem 3.2: Let ^ be a simple graph with Laplacian 
matrix Lg. Then, the first three spectral moments of the 
Laplacian matrix are 

TOi {Lg) = -Si, (3) 

n 

"^2 {Lg) = - {Si + 5*2) , 
n 

"^3 (is) = -(3^2 + 53 -6 A), 
n 

where Sp is defined in (fTh, and A is the total number of 

triangles in Q. 



Proof: First, we apply Lemma 3.1 to compute the Lapla 



cian spectral moments of a simple graph Q in terms of closed 
walks in the weighted Laplacian graph C {Q), described in 



Definition 3.1 In particular, the first spectral moment is related 
to the set of closed walks of length 1. Clearly, the only possible 
closed walks of length 1 in £ {Q) are those using a self-loop, 
i.e., Pl]]^ — {{vi,Vi) for Vi e V}. Since the weight uj {p) of 
a walk p ~ {vi, vi) in C {Q) is equal to di, the first spectral 
moment is equal to: 



nil {Lg) 



1 



viev 



Uj{{Vi,Vi)) 



n 

n ^-^ 

4 = 1 



1 



-Si. 



In order to compute the second moment ni2 {Lg), we need 
to account for two different types of closed walks of length 
two: 

(a) Walks in C {G) of the first type visit a self-loop twice, 
i.e., p = {vi,Vi,Vi), and the weight of this walk is d^. We 
denote the set containing walks of this type as P2J . 

(b) Walks of the second type start at a node Vi, visit a 
neighbor Vj e Afi, and return to Vi. This type of walk defines 



the set Pj^ = {{vi, Vj, Vi) for Vi £ V, Vj G A/i}. The weight 
of any walk of this type is equal to 1. 

Since p!^ U P^ = p!^l and P^} n p!^ = 0, the sets p!^ 



,(0 



3(*) 



2 „ aiiu J 2a ' ' ^ 2b 

and P^'f^ form a partition of the set of walks P2n- Hence, 
from Lemma 13. 11 we have 



m2{Lc, 



n ^ ^ 

1 V-.2 



■b) 



1 



E E' 



■b) 



ri ^ ^ ^ ^ 



Vi&> 



vi^V vjeMi 



1^-^,2^! Vd, 



i=l 






- (52 + Si) 
n 



In order to compute the third Laplacian spectral moment, 
we classify closed walks of length 3 in £ (Q) into sets that 
partition P3 „. In particular, we classify these closed walks 
with respect to the structure of the subgraph covered by the 
walk. Specifically, two walks pi and p2 belong to the same 
type if the subgraphs covered by the w alks, denoted by C (pi) 
and C {P2) according to Definition 3.2 are isomorphicn 
In Fig. |2] we depict the set of all possible nonisomorphic 
subgraphs covered by closed walks of length 3 (in solid lines). 
According to this classification, we have the following types 
of walks: 

(a) Walks of the first type present the form {vi,Vi,Vi,Vi), 
covering the self-loop three times. The subgraph covered 
by this walk is the isolated node, represented in Fig. l2|a). 
Although the walk covering this graph visits the self-loop 
{vi, Vi} (in dashed line in Fig. l2fa)), this self-loop is not part 
of the covered subgraph (according to Definition |3.2|i. The 



,Vi) 



weights associated with walks of this type are equal to df 
and we denote the set of this type of walks as P^J . 

(b) Walks of the second type present the form {vi , v 
and {vi,Vj,Vj,Vi) with Vj E Ni- The subgraph covered by 
these walks is the single edge in Fig. l2|b), where the self- 
loops (in dashed line) are used by the walks but not included 
in the subgraph. These two walks have weights di and dj, 
respectively, and are represented in Fig. |2|b). We denote this 



set of walks as P. 



(*) 



36 



(c) The last type of walk presents the form {vi,Vj,Vk,Vi) 



and {vi,Vk,Vj,Vi) for v^ 



Vk 



Vi. The subgraph 



covered by these walks is the triangle (Fig. l2|c)). These walks 
have weights —1 and we denote this set of walks as P, 



j{i) p(«) 



Since the sets P^J, P3V' ^"d P^J partition the set of closed 



(») 



(0 

3c ■ 



(i) 

walks P3 „, we have that 



maiLg) 



E E 



E ^(py 



v,ev xe{a,b,c} p^pg) 



We can now analyze each of the terms in the above summation 

For convenience, we define T^r^ — ^ 

and study each one of the terms T^^ for x G {a,b, c 



^Et,,evEpgp(o w(p)' 



'*Two graphs G and H are isomorphic if there is a bijection between the 
vertex sets V{G) and V{H), f : V{G) — )■ V(H), such that any two vertices 
u and V oi G are adjacenct in G if and only if f{u) and f{v) are adjacent 
in H. 



4.,--. 



(a) 



(b) 



"j 




Fig. 2. Collection of possible subgraphs covered by closed walks of length 3. 
Since these subgraphs are simple, by definition, they do not contain self-loops 
(which are included in the figure for convenience in our derivations). 



(i) 

(a) For X = a, all the walks in Pg^ have a weight equal to 
dl Hence, we have Tg^ - i ^7=1^1 = Ss/n. 

(b) For X — b, the walks in Pg^ have a weight equal to 
either d, or do. Then, we have: 



T3b 



1 



T. d^ + ^T. E dj 



n n 



dij Cli 



1 A .0 1 



^ n n 

- E E "'j^j 



i=i 3=1 



n 



E^^ 



Y^ j2 



df + -2^d'j= 2S2/n. 



i=l 



n 



j=i 



(c) For a; = c, we have Ta^ = ^ Ei<»,i,fe<„ ("1) *yfe' 
where i,,i. is an indicator function that takes the value 1 



when Vi 



''ijk 



Vk 



Vi (0, otherwise). Hence, Tj 



3c 



n Si<i<7<fc<n (~1) ^yfc' where the coefficient 6 accounts for 
all the possible permutations of the three indexes i,j, and k. 
Hence, since J2i<i<i<k<n^ijk = A/n (where A is the total 
number of triangles), we have that T^c — — 6A/?7. 

Finally, since 7713 {Lg) — T^a + T-^b + P3C, we obtain 
the expression for the third moment in the statement of the 
theorem. ■ 



Remark 3.2: Theorem 3.2 relates purely algebraic proper- 
ties - the spectral moments - to structural features of the 
network, namely the degree sequence and the number of 
triangles. The key concept behind the proof is the classification 
of closed walks into sets according the subgraph covered 
by these walks. This idea can be extended to higher-order 
moments, although the combinatorial analysis becomes more 
elaborate as we increase the order of the moments. 

Remark 3.3: Notice that the power-sums of the degrees, 
Sp, and the number of triangles. A, can be retrieved from 
the set of first-order neighborhoods, {Ql for all Vi G V}. 
Therefore, we can compute the first three Laplacian spectral 
moments from local structural measurements with radius I. 
As we shall discuss below, it is possible to compute the first 
2r + 1 Laplacian spectral moments of a graph whenever we 
have access to the set of r-th order neighborhoods, {Qi for 
all v,GV}. 

2) Higher-Order Laplacian Spectral Moments: We now 
extend our analysis to the fourth- and fifth-order spectral 
moments of the Laplacian matrix. We first define the col- 
lection of structural measurements that are involved in our 
expressions. Let us denote by ti, qi, and pi the number of 
triangles, quadrangles, and pentagons touching node Vi in Q, 
respectively. The total number of quadrangles and pentagons in 



Q are denoted by Q and P, respectively. The following terms 
define structural correlations that are relevant in our analysis: 



f-^dd — 



-Y^dd,, Cd2d = -y2 d^.dj, (4) 

77 ^ ^ 77 -^ ^ 



Cdt = — y ^ diU, Cd^t = — y d^ti, 

77 ^ ^ 77 ^ ^ 



Ti.ev 



w.ev 



Orf„ — 



- V d,q,, Ddd^- V d,dj |M n AAjI , 

77 ^ ^ 77 -^ ^ 
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where |A/i n JVj \ is the number of common neighbors shared 
by Vi and Vj. The above terms represent correlations between 
structural variables in the graph, some of which have already 
been studied extensively. For example, Cdd and Cd^d are 
closely related to the widely studied assortativity coefficient 
lis], which measures the preference (or aversion) of nodes to 
attach to nodes presenting a similar degree. On the other hand, 
we find in (|4]) several correlation terms that have not been 
studied in such depth. For example, the terms Cdt and Cd^t 
quantify the correlation between degrees and local clustering^ 
throughout the set of nodes. Although nontrivial variations of 
local clustering with respect to degrees have been reported 
in the complex networks literature ifTSl . the effects in the 
network's behavior are not well understood. 

The main result in this section relates the fourth and fifth 
Laplacian spectral moments to local structural measurements 
and correlation terms, as follows: 

Theorem 3.3: Let Q he a simple graph with Laplacian 
matrix Lg . Then, the fourth and fifth Laplacian moments can 
be written as 



fni{Lg) = 



1715 (Lg) 



n 
+^Cdd - 

- (-5^2 



-10 {C, 



dd 



- 25*2 + 453 + S4 
&Cdt, 
+ 55-34 
+ Cd'^d 



5544 

— Cdt 



c. 



d^t 



(5) 



30A - lOP) 

Cdq — Ddd) 



where Sp is defined in (fill, and the correlation terms Cdd, Cdt, 
Cdq, Cd2d, Cd2t, and Ddd are defined in (HI). 



Proof: As in Theorem 3.2 the proof is based on a 
classification of closed walks of length 4 and 5 with respect 
to the subgraph covered by the walks. Details about this 
classification and the subsequent combinatorial analysis can 
be found in the Appendix. ■ 

Remark 3.4: Observe how, as we increase the order of the 
moments, more complicated structural features arise in our 
expressions. In particular, the fifth moment is influenced by 
the degree sequence, the number of cycles of length 3 and 5, 
and all the correlation terms defined in (|4]). 

Remark 3.5: As we mentioned before, the maximum order 
of the spectral moment that one can compute from local struc- 
tural information depends on the radius r of the neighborhoods 
that are accessible to us. One can prove that in order to 
compute the fc-th Laplacian spectral moment, we need to count 

^The local clu.stering coefficient is defined as t; / ( ') , (5]. 



the number of cycles of length k in Q. One can also prove 
that access to neighborhoods of radius r allows us to count all 
cycles of length up to 2r+l in Q, fT^l. Hence, we can compute 
the set of Laplacian moments up to the order fcmax = 27- + 1 
using local structural information that can be extracted from 
the set of r-th order neighborhoods, {Gi for all Vi e V}. 

In this section, we have derived expressions to compute 
the first five spectral moment of the Laplacian matrix of 
a network from local structural measurements. In the next 
section, we present a series of semidefinite programs (SDP's) 
whose solutions provide optimal bounds on the Laplacian 
spectral radius and spectral gap. 

IV. Optimal Laplacian Bounds from Spectral 
Moments 

In this section, we introduce a novel approach to compute 
bounds on the spectral gap and the spectral radius of the 
Laplacian matrix from a truncated sequence of Laplacian 
spectral moments. More explicitly, the problem solved in this 
section can be stated as follows: 

Problem 1 (Moment-based bounds): Given a truncated se- 
quence of Laplacian spectral moments (m^ {Lg))^^^, find 
bounds on the spectral gap and the spectral radius of the 
Laplacian matrix Lg. 

Our results are based on an optimization framework recently 
proposed in 1*201. In order to adapt our problem to this 
framework, we need to introduce some definitions. Given a 
simple connected graph Q with Laplacian eigenvalues {Ai}"^j^, 
we define the spectral density of the nontrivial eigenvalue 
spectrum as 



Pe(A) 



^^<5(A-A,), 



(6) 



where 5{-) is the Dirac delta function. Notice how we have 
excluded the trivial eigenvalue, Ai := 0, from the spectral 
density; hence, the supportFlof pg (A) is equal to supp{pg) ~ 
{Ai}"^2- The moments of the spectral density in (raj), denoted 
by mfc (Lg), can be written in terms of the spectral moments 
of Lg, as follows 

lukiLg)^ / A'^ -V5(A-A,) d\ 



i=2 



x' = 



n-\ 



■mk(Lg), 



(7) 



for all fc > 1 (where we have used the fact that Ai = 0, in 
our derivations). 

In what follows, we propose a solution to Problem [T] using 
a technique proposed by Lasserre in ll20l . In that paper, the 
following problem was addressed: 

Problem 2: Consider a truncated sequence of moments 
{^Ir)i<:r<k corresponding to an unknown density function 
/i(A). Denote by [a, 6] the smallest interval containing the 

^Recall that the support of a density function ^ on R, denoted by supp (fi), 
is the smallest closed set B such that ^ (R\_B) = 0. 



support of fi. Compute an upper bound a > a and a lower 
bound /3 < 6 when only the truncated sequence of moments 
is available. 

In the context of our spectral problem, we have a truncated 
sequence of five spectral moments, (mr(ig))i<r<5, corre- 
sponding to the unknown density function pg and given by 
the expressions ([3]l, dSll, and (JTl). In this context, the smallest 
interval [a, b] containing supp{pg) is equal to [A2, A„]. There- 
fore, a solution to Problem [2] would directly provide an upper 
bound on the spectral gap, a > X2, and a lower bound on 
the spectral radius, /3 < A„. We now describe a numerical 
scheme proposed by Lasserre in fSOl to solve Problem [2] This 
solution is based on a series of semidefinite programs in one 
variable. In order to formulate this series of SDP's, we need 
to introduce some definitions. For any s e N, let us consider 
a truncated sequence of moments M = iMr)/^^ , associated 
with an unknown density function p. We define the following 
Hankel matrices of moments: 



i?2. (M) ^ 



R 



2s+l 



(M)^ 



1 


Ml •• 


• M, 


Ml 


M2 •• 


• M,+i 


M. 


M.+i •• 


■ M2s 


Ml 


M2 


■ ■ ■ Ms+i 


M2 


M3 


■ ■ ■ Ms+2 



M,+i M,+2 • • • Af2,+i _ 
We also define the localizing matriiP\Hs (a;,M) as, 
H, {x, M) ^ i?2s+i (M) - X R2s (M) . 



(8) 



(9) 



(10) 



The localizing matrix presents the following two properties 
ED: 

PI. Hs {a, M) >: 0, if the support of fi is contained in 

the set [a, 00), and 
P2. —Hs {(3, M) ^ 0, if the support of fi is contained in 

the set (— C)o,/3]. 
Based on these properties, Lasserre proposed the following 
series of SDP's to find a solution for Problem |2] 

Solution to Problem 2: Let M = (-^•^r)rli bs a truncated 
sequence of moments associated with an unknown density 
function p. Then 

a < a, (M) = max {a; : iJ^ (a;, M) ^ 0} , (11) 

X 

b> l3s{M)^mm{x:~Hsix,M)'^0}, (12) 

X 

where [a, b] is the smallest interval containing the support of 

p. 

Remark 4.1: Note that the entries of the localizing matrix 
Hs {x, M) depend affinely on the decision variable x. Thus, 
as (M) and /3s (M) are the solutions to two SDP's in one 
variable, which can be efficiently solved using standard opti- 
mization software. 

'a more general definition of localizing matrix can be found in 121]. For 
simplicity, we restrict our definition to the particular form used in our problem. 



Therefore, we can directly apply the above result to solve 
Problem [ll by considering the sequence of moments m = 
{mr{Lg))/J'^ = {i^mr {Lg))l''Jj^ in the statement of 
the solution to Problem l2] Since this sequence of moments 
corresponds to the spectral density pg, with support {Ai}"^2' 



the solutions in (111 and (12 1 directly provide the following 
bounds on the spectral radius and spectral gap: 

Solution to Problem 1: Let m = {^p^rrir (ie))r=V ^^ 
a truncated sequence of (scaled) Laplacian spectral moments 
associated with a graph Q. Then the Laplacian spectral gap 
and spectral radius of Q satisfy the following bounds: 

A2 < cis (Hi) = max {a; : Hs {x,m) >: 0} , (13) 

X 

A„ > f3s (m) ^ min {x : -Hs {x, m) ^ 0} . (14) 



In Section III-B| we derived expressions for the first five 
Laplacian spectral moments, {nir (ie))r=i> in terms of local 
structural measurements of the network. Therefore, we can 
apply the above solution to find spectral bounds for s = 2. 
Furthermore, for s = 2 we can prove that the optimal values 
a2 (in) and /32 (iii) are the maximum and minimum roots of a 
cubic polynomial, as follows. For s = 2, the localizing matri?!^ 
takes the form 





mi — X 


7112 


— xWli 


rn^ — xfn2 


2 {x, m) = 


rfi2 — xmi 


ma 


— xfn2 


Wii — xfn^ 




7713 — xfiT,2 


m4 


— xfriT, 


Trie, — a;m4 



Let us first analyze the optimal value /32 (iii) in (14i. Note 
that — _ff2(a;,iii) ^ if and only if all the eigenvalues of 
H2 {x, m) are nonpositive. The characteristic polynomial of 
H2 {x, in) can be written as, 

02(A) ^ det(A/-i?2(a;,rii)) 

= A^ + ri (x) A^ + r2 (x) X + r^ (x) , 

where rj (x) is a polynomial of degree j in the variable 
X (with coefficients depending on the sequence of moments 
{Wii, ...,rn^)). One can prove that all the eigenvalues of the 
Hankel matrix H2 (a;,iii) are real ||22I . Also, by Descartes' 
rule, all the roots of 02 (A) are nonpositive if and only if 
T^j {x) ^ 0, for j = 1,2, and 3. Lasserre proved in EOl 
that rj (x) > if and only if r^ (x) > 0. Since r^ (x) = 
— detH2 (a;,m), we have that the optimal value /32 (rri) in 



(14 1 is equal to the largest root of the cubic polynomial 
ra (x) = — dctiJ2 (a^jiii) — 0. Similarly, one can also prove 



that the optimal value a2 (m) in (13i is the smallest root of 
ra (x). We now state this result more explicitly as follows: 

Solution to Problem 1 with 5 Moments: Let m = 
i^ri^^^^ i^o))r=i ^^ ^ truncated sequence of (scaled) Lapla- 
cian spectral moments associated with a graph Q. Then the 
Laplacian spectral gap and spectral radius of Q satisfy 

A2 < a2 (in) = min{xi,X2,a;3}, (15) 

A„ > ^2 (iii) =max{xi,a;2,a::3}, 

where a;i,a;2, and 2:3 are the roots of ra (x) = 

— det i?2 (a;, m). 

*We shall remove the argument Lg from the Laplacian spectral moments 
Trir {Lg), for simplicity in notation. 




/, cycle length 
(c) 

Fig. 3. In (a), we plot the complementai'y cumulative distribution of degrees for the electrical transmission network under study. In (b), we plot the joint-degree 
distribution J^j (fci, ^2) in the range I < ki,k2 < 9 for the American power grid. In (c), we plot on a semilogarithmic scale the number of cycles (pi in the 
thi'ee transmission networks for lengths 3 < I < 10. 



Remark 4.2: Note that it is easy, but tedious, to write down 
an explicit expression for r^ (x) in terms of the Laplacian 
spectral moments. Also, we can find explicit expressions for 
the roots of the cubic polynomial r^ (x) using, for example, 
Cardano's rule t23]. Therefore, it is possible to derive explicit 
expressions for 02 (m) and 13-2 (m) in terms of the Laplacian 
spectral moments, but the resulting expressions are so com- 
plicated that they do not provide much insight. 

In this section, we have presented an optimization-based 
approach to compute optimal bounds on the Laplacian spec- 
tral gap and spectral radius from a truncated sequence of 
Laplacian spectral moments. These spectral moments can be 
written in terms of local structural measurements using ([3]) 
and dSll. Hence, our results build a bridge between structural 
measurements of a network and its spectral properties. In the 
following section, we use our theoretical results to analyze 
spectral properties of electrical transmission networks. 

V. Structural Analysis and Simulations 

Motivated by recent interest for smart grid architectures, 
there is a fast-growing literature studying the structure and 
function of electrical transmission and distribution networks 
ll24l - ||26ll . Most of the results found in the literature make 
use of extensive numerical simulations to find relationships 
between structural properties of a power grid and the behavior 
of dynamical processes taking place within it. In this section, 
we apply our theoretical results to study the relationship 
between structural and spectral properties of the American 
(western states), French and Spanish transmission networks. 
Our analysis reveals that the Laplacian spectral radii and spec- 
tral moments of electrical transmission networks are strongly 
constrained by the degree sequence and the so-called joint- 
degree distribution. We also verify that the spectral gap cannot 
be efficiently bounded using local structural features only, 
since this spectral property strongly depends on the global 
connectivity of the network. 

A. Local Structural Analysis 

We examine unweighted, undirected graphs representing the 
structure of high-voltage transmission networks in the USA 



(western states), France and Spain (the adjacencies of these 
networks are available, in MATLAB format, in p27l|). The 
number of nodes (buses) and edges (transmission lines) in 
these networks are: n^^ — 4, 941 and e^g — 6, 594 for the 
American network, 71/,. = 146 and ejr = 223 for the French, 
and Tisp = 98 and e^p = 175 for the Spanish. From these data 
sets, we compute a collection of structural features that are 
relevant to the network's functionality. 

We begin our structural analysis by studying the distribution 
of degrees in the network. In the network science literature, 
it is common to classify the degree distribution by looking at 
the complementary cumulative distribution function (CCDF) 
of the degrees, denoted by ccdf(fc). The function ccdf(fc) is 
defined as the probability that a node chosen uniformly at 
random presents a degree higher than or equal to a given 
value k. In Fig. [3ta) we plot the CCDF's for the three 
transmission networks under study. These distributions are 
well approximated by exponential functions |24|. 

Another extensively studied structural property is the joint- 
degree distribution. This distribution is related to the number 
of edges connecting sets of nodes of different degrees |28l, 
and it is defined as follows: Consider the sets of nodes with 
degrees fci and k2, denoted by Vk^ and Vfcj; then, the joint- 
degree distribution is defined as: 

where e is the total number of edges in the network. The 
correlation terms Cdd and Cd^d defined in Q can be explicitly 
related to the joint-degree distribution, as follows 



C, 



dd 



- V" V" kik2Jd{ki,k2) , 
n ^-^ ^-^ 

ki = l fc2 = l 
dmax dmax 

Cd^d = -"^ "^ klk2Jd{ki,k2), 



ki = l fc2 = l 



where dmax is the network's maximum degree. In Fig.lSlb), we 
plot a contour plot representing the joint-degree distribution of 
the American power grid (for 1 < ki,k2 < 9). 



TABLE I 

SUMMANDS INVOLVED IN THE FIRST FIVE LAPLACIAN MOMENTS OF THE AMERICAN GRID. 
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- 
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- 
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742.4 


-2.66 


20.66 


229.1 


432.4 


- 


- 


1.58 


- 


80.12 


- 


-18.81 


- 


- 


- 


5 


7,588 


- 


-51.6 


286.4 


2,162 


4,157 


3.95 


- 


-3.68 


200.3 


1,127 


-23.51 


-177.2 


5.30 


14.58 



A structural distribution that is relevant in our analysis is 
the number of cycles in the network (i.e., closed walks with no 
repeated nodes). We denote by (/>/ the total number of cycles 
of length I in the network, and plot (/>/ versus I for the three 
networks under consideration in Fig. l3jc). 

In previous sections, we have provided tools that allow 
us to quantify the effects of these (and other) structural 
measurements on the Laplacian spectrum. We base our subse- 
quent analysis on the closed-form expressions for the spectral 
moments in (|3]l and (|5]l. Note that these expressions can be 
written in a unified form as follows: 



rrik {Lq 



t^ 



(r) ^r 

n 



s=3 



n 



t=i 



(16) 



where Sr is the r-th power sum of the degrees, (fig is the 
number of cycles of length s, and Ct indicate the correlation 
iCdd,Cd^dTCdt,Cd2t,Cdq,Dj_ci)- The 
are, respectively, the coeffi- 



terms with (C()j 
coefficients S 



(r) 






and %' 



cients accompanying Sr/n, tps/n, and Ct in the expression 
for the k-th moment. For example, according to p), the third 

;3r_ 



spectral moment has coefficients S. 



(2) _ 



3, S. 



1, and 



73 



(3) 



-6 (and the rest of coefficients are zero). 



We now use ( 16 1 to quantify the effect of relevant structural 



measurements on the spectral moments of the American 
power grid. In this network, the relevant structural measure- 
ments are: (/) the power-sums of the degrees: {Sr)r^i ~ 
(2.669,10.33,57.29,432.4,4157), (ii) the number of cycles: 
03 = A = 651, 04 = Q = 979, (f)^ = P = 1821, and (Hi) the 
correlation terms Cdd = 20.03, Cd^d = 112.7, Cdt = 2.35, 
Cd2t = 17.72, Cdq = 5.30, and Ddd = 14.58. In TableJ!] 
we include the values of each one of the summands in ( |16^ 
for the first five Laplacian moments of the American power 
grid. In this table, we can detect those terms that dominate 
the sum for each spectral moment. We have marked (in bold 
numerals) those summands that individually account for over 
10% of the total sum; the sum of these dominant terms in each 
case provides a very good approximation to nik- We observe 
how the power-sums of the degree sequence Sr have a large 
impact on all the spectral moments. Also, the correlation terms 
Ci = Cdd and C2 = Cd^d are significantly important for the 
fourth and fifth moments. 

From these empirical observations, we conclude that the 
degree sequence plays a key role in the Laplacian spectral 
moments, as expected. Furthermore, the correlation terms Cdd 
and Cd^d also have an important influence on the spectral 
moments. Note that these correlation terms, defined in Q, 
directly depend on the joint-degree distribution, Jd{ki,k2)- 
We conclude that the structural measurements that are most 
relevant for the Laplacian eigenvalue spectrum are: 



(/) The degree distribution, via the power-sums Sr- 
(ii) The joint-degree distribution, via the correlation terms 
Cdd and Cd^d- 

Other correlation terms and the distribution of cycles of 
lengths 3 to 5 have a weak influence on the spectral moments. 
Similar structural results are also observed for both the French 
and Spanish transmission networks. The methodology intro- 
duced here allows us to recognize the degree distribution and 
the joint-degree distribution as the structural properties that 
are most influential on the Laplacian spectraL moments of the 
electrical grids under study. 

B. Spectral Bounds from Structural Information 

From the collection of structural measurements studied 
above, we compute the first five spectral moments for the 
transmission networks under study, seeTable lU] We then use 
these moments to compute bounds on the spectral gap and 
spectral radius, a-2 and P2, using the methodology described 
in Solution to Problem 1 with 5 Moments, in Section |IV] 
The numerical values for these bounds, as well as the exact 
values for the spectral gap and spectral radius, A2 and A„, 
are included in Table [III] In this table, we also include the 
values obtained for the bounds in ||T3| and |fT6l (which were 



described 


in Section III-A|i. 

TABLE II 
First five Laplacian moments. 
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m2 (La) 


"13 [Lg) 
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ms (Lg) 


USA 
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Spain 


2.66 
3.05 

3.57 


13.00 
15.58 
20.83 


87.50 
97.73 
147.3 


742.4 
680.1 
L155 


7,588 
5,046 
9,686 



TABLE III 

Upper and lower bounds on the Laplacian spectral gap and 
spectral radius, a2 and a„. 



Network 


113 


A2 (Lg) 


02 


^2 


An (Lg) 


l'161 


USA 
France 
Spain 


1.16e-4 
3.70e-3 
4.90e-3 


7.59e-4 
4.24e-2 
7.70e-2 


0.97 
0.81 
0.86 


12.8 
8.07 
9.18 


20.10 
9.60 
10.66 


22.28 
9.86 
10.79 



In Table [III] we observe that /32 and the bound in lfT6l 
strongly constrain the intervals of possible values for the 
spectral radii, specially in the case of European transmis- 
sion networks. In Fig. |4|a)-(c), we plot the histograms for 
the Laplacian eigenvalues of the American (western states), 
French and Spanish transmission networks, where we have 
also included the bounds in Table [III] We observe that 
bounds for the Laplacian spectral gap based on local structural 
measurements perform poorly in practice. The main reason 



10 



us Power Grid 

3r ^ 




French Power Grid 



(X2=0.97 



3 




?.i=0 X2=4.2e-2 

i 1 


y^ 0.02 0.04 0.06 0.08 

j=3.7e-3 

X„=9.60 


^^^^^" 


6 ;r t 



|32=8.07 p=9.86 



(b) 



Spanish Power Grid 




8 Tio f 

P2=9.18 p=10.79 



Fig. 4. In Figs, (a), (b) and (c), we plot liistograms for the Laplacian eigenvalues of the American (western states), French and Spanish liigh-voltage 
transmission networks, respectively. In the small subfigures, we plot a histograms of the eigenvalues around the origin and indicate the spectral gap. We also 
include the values of the bounds of the Laplacian spectral radius and gap from Table [ill] 



behind this Umitation is that the Laplacian spectral gap is a 
global property that quantifies how 'well-connected' a network 
is Il29l .ll30l. Since the structural measurements used in our 
bounds (degree sequence, correlation terms, etc.) have a local 
nature, they do not contain enough information to determine 
how well connected the network is globally. 

C. Spectrum-Preserving Structural Perturbations 

In this subsection, we validate our previous structural analy- 
sis with numerical computations. One of our main conclusions 
is that the degree and joint-degree distributions are the struc- 
tural features with the strongest influence on the Laplacian 
spectral moments and spectral radii of electrical transmis- 
sion networks. We can empirically verify this conclusion by 
studying the effect of structural perturbations that preserve 
these features on the Laplacian spectrum. In what follows, we 
provide an algorithmic description of a structural perturbation 
that transforms an input graph Gin — (V,£in) into a different 
output graph Gout = (V,£out) while preserving both the 
degree and the joint-degree distributions; 

(/) Randomly choose two degrees fci and ^2 drawn from 
the joint-degree probability distribution Jd {ki, k2). 

(ii) Choose two edges {vi,Vj} and {vr,Vs} uniformly at 
random from the set of edges in £in that satisfy d,; = d,. = fci 
and dj — ds — k2. 

(Hi) Build an intermediate network F — (V,£+), where 

5+ = Sin + {Vi,Vs} + {Vr,Vj} - {Vi,Vj} - {Vr,Vs} , 

i.e., remove edges {vi,Vj} and {vr,Vs}, and add edges 
{vi,Vs} and {vr,Vj}. 

(iv) If F is connected and has no multi-edges (several edges 
connecting the same pair of nodes), we define the output 
rewired network Gout — F. Else, go to step (/). 

Note that the above algorithm defines a random rewiring 
that maintains both the degree distribution and the joint- 
degree distribution of Gin, since the number edges between 
the sets of nodes with degrees fci and ^2, denoted by Vk^ and 
Vfcj, is not modified. Since these structural properties have 
the strongest influence on the Laplacian spectral moments, 
we should expect the Laplacian spectrum to remain almost 
invariant. In order to verify this statement, we recursively 
apply this random rewiring 5, 000 times to the American 



power grid. We can verify that, despite the large number 
of rewirings, the spectral radii and spectral moments of the 
resulting network, which we denote by Qrw, are practically 
identical to those of the American power grid. In particular, 
the percentage of change in the first five spectral moments of 
Qrw is under 0.5%. Furthermore, the spectral radius of Qrw 
is 20.18, which represents a 0.36% change with respect to 
the original network. As expected, the spectral property that 
is most sensitive to the rewirings is the spectral gap, which 
takes a resulting value of 0.016 (from a value of A2 = 8.4e — 4 
for the American grid). 

VI. Conclusions 

This paper studies the relationship between local structural 
features of large complex networks and global spectral prop- 
erties of their Laplacian matrices. We have proposed a novel 
approach, based on algebraic graph theory and convex opti- 
mization, that allows us to quantify the effect of an important 
collection of local structural features on the Laplacian spectral 
moments and spectral radius. The following is a list of our 
main contributions and conclusions: 



1) We have derived in Section pII-B| explicit expressions 
to compute the first 2k + 1 spectral moments of the 
Laplacian matrix using local structural measurements 
with radius r. 

2) We have proposed an optimization framework in Section 
[TV] to compute lower bounds on the spectral radius of 
the Laplacian matrix of a network from a truncated 
sequence of spectral moments. Our bounds take into 
account the effect of important structural properties that 
are usually neglected in most of the bounds found in the 
literature, such as the distribution of cycles and structural 
correlations. 

3) Our analysis reveals that the local structural measure- 
ments that are most relevant to the Laplacian spectrum of 
electrical transmission networks are the degree sequence 
and joint-degree distributions. We have numerically 
verified that other structural properties that may seem 
important, such as the distribution of short cycles in the 
network, have a very weak influence on the Laplacian 
spectrum. 
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4) As expected, local structural features are not enough 
to efficiently bound or estimate the Laplacian spectral 
gap, since this spectral property strongly depends on the 
global connectivity of the network. Since the spectral 
gap is fundamental in the analysis of many dynamical 
processes on networks, random models in which only 
local structural features are prescribed, such as those in 
Il7l-f9l, are not valid to generate synthetic topologies in 
which these dynamical processes can be studied. 
Although the mathematical tools developed here can be ex- 
tended to the case of weighted networks, we have limited our 
analysis to unweighted networks to isolate the role played by 
the network topological structure from other factors. Ongoing 
work addresses the spectral analysis of large-scale networks 
with weighted edges, as well as networks with directed edges. 
Finally, the analytical approach introduced is also applicable 
to the analysis of spectral properties of many real networks in 
other areas of science and technology, such as networks of 
multi-agent systems, communication networks, as well as so- 
cial and biological networks. Since the Laplacian eigenvalues 
are relevant in a wide variety of dynamical processes , our 
approach is useful to theoretically analyze the effect of local 
structural properties on those processes. 

Acknowledgement 1: We thank V. Rosato for kindly pro- 
viding us with the topologies of the French and Spanish 
transmission networks studied in 



Appendix 



Theorem 3.3 Let tj be a simple graph with Laplacian matrix 
Lg. Then, the fourth and fifth Laplacian moments can be 
written as 



TO4 {Lg) 



mz{Lg) = 



1 



{-Si + 25*2 + 453 + 5*4 

-4:Cdd — &Cdt, 



1 



- (-5^2 - 
n 

+10 {Cm 



55.3 + 554 + 55 



Cd'^d ^ Cdt — Cd^t 



30A - lOP) 

Cdq — Ddd) 

where Sr = X^u ev "^I' ^^^ '^he correlation terms C^d, Cdt, 
Cdq, Cd2d, Cd2t, and Ddd are defined in Q. 

Proof: As in Theorem 3.2 we use Lemma [TT] to compute 
the Laplacian spectral moments in terms of weighted sums of 
closed walks in the weighted Laplacian graph Cq. In order 
to compute the fourth Laplacian spectral moment, we classify 
the types of possible closed walks of length 4 into subsets 
according to the structure of the underlying graph covered by 
the walk. As we explained in the proof of Theorem 3.2 two 



walks of length 4 are in the same subset if the subgraphs 

they cover are isomorphic (without considering self-loops). 

We enumerate the possible types in Fig. Island we denote the 

corresponding sets of walks as P^l^, P^^', P^l', P^'J, and PJ^ . 

(i) 



,(«) 



These sets P\^J^,...,PlJ partition the set of closed walks P^^ 
Hence, we have 

■mi{Lg) 



Vi^V a;G{a,b,c,d,e} p^p'^^) 



We now analyze each one of 
above summations. For convenience. 



the 



terms in 
define T^x 



the 



!"'\ 


1 \ 


-' 






V/ 




y 


v 


-1 \/ 


^k* 


y 


/^ 


< 


V,- 


'-. 


^v. 



id) 



ib) 



(c) 



d< ,— . 



dk ,— - 




-1 



.4 



( !^k -1 Vj\ ) V, 



-1 



Fig. 5. Collection of possible graphs covered by closed walks of length 4. 

^Lt, evE„pp(')^(p) and analyze the term T^^ for x € 
{a, b, c,d}: 

(a) For a; = a, we have that the weights uj {p) of the walks 
in P|*^ are all the same, and equal to df. Hence 



T4<j = - V'df = 54/r 

77 -^ ^ 



(b) For X — b, the weights of the walks in P|^^ are equal 

to 2 + 4 (df -I- d| + didj). Hence 



T4b 






2 + 4 (d- + 4 + didj) 



1 



(5i+453)+4C, 



dd- 



(0 



(c) For X = c, the weights of the walks in P^J (i.e., walks 
that cover the two-chain graph) are equal to 4. Hence 



n, 



I E ^'■2itfl'V = §(s.-5. 



n ^ — ' n ^ — ' \ 2 



where in equality (/) we have used the fact that the number 

of two-chain graphs whose center node is Vi is equal to (2'). 

id) For X ^ d, the weights of the walks in P^^ are equal 

to —8 {di + dj + dk). Hence, 



T4,d 



Y, -8{d,+dj+dk) 



n n n 



8 



^^EEE^^^j'^^^' 



i=i j=i k=i 

where tijk is an indicator function that takes value \ \f Vi ^ 
Vj ^ Vk ^ Vi. Since YJ^=i YJk=i ^^ijk = U (the number of 
triangles touching node Vi), we have that 



lid — / , i-i"! 

n 

2 = 1 



Tid — ^_^tidi — —8Cdt- 



,(«) 



(e) For X = e, the weights of the walks in P^J are equal 
to 8. Hence, 

1 



T4, 



E ' 

s.t. l<i<j<A:<r<n 



8Q/n. 
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Finally, since TO4 (Lg) = T4a +T4b + T4c + T4d, we obtain where 6^ = X]fc=i '^jfe'^ife = \J\fi O J\fj\, the number of 
the expression for the fourth Laplacian spectral moment in the common neighbors shared by Vi and Vj . Hence, 
statement of the theorem after simple algebraic simplifications. 

In order to derive a similar expression for the fifth-order 30A 10 ^^ 72 10 V^ 

Laplacian spectral moments, we follow an identical approach. 5d — ^ ^ 2_^ * « ^ Z_^ *J «J \ ^ j) 

In this case, the algebraic manipulations become more tedious. *~ *~"' 

Below, we provide the main steps in the derivations (and omit ^ '" ''^* '*'' 

the details regarding the algebraic manipulations). As before, 

we partition the set of closed walks P^l according to the ^^) Fo^ ^ ^ ^^ *e weights of walks covering the quadrangle 

subgraph covered by the walk. We show the structure of the S^aph are d,, d,, d^, and d^. Countmg the multiphcities of each 
possible subgraphs in Fig. |6] 'yP^ °f ^^^^ we have that 

We now analyze each one of the terms T^^ — 1 

^E.,evE,ep'')^(p)fo'-xe{a,5,...,g}: ^^e = - Y. W (d, + d, + d^ + dr) 

(a) For X = a, we have 



n n n n 



■i— 1 J — 1 fc— 1 r— 1 



T5a = -E^'=^5Ai. 



n 

i=l 



... „ , , . „ ., , , , ,, where Qiikr is an mdicator function that takes value 1 if 

(b) For X = 0, we can determine all possible closed walks ■' r,- ■s-^ ^ 

f 1 fu 5 • .u ^ u • c- [zvu\ ^ ^ • .u . i;^ ~ ■Uj -- Wfe ^ Wr ^ Wj- Since Z^k,- j. r<n4gij'cr = 9* 

of length 5 using the edge graph in Fig. 6[b) and derive that , , ■' , j. , , , . ^^J^i:J^'^'''i:" / 

■— 1^ (the number of quadrangles touching node Vi), we have that 



T56 = - E ^ (^^^ + '^1 + rff + 4 + c^* cij + dtd]) ^Q n 

Vi~Vj T^e — > qidi = lOCdg. 

= - (52 + 54) + 10Cd2rf. 

(/■) For a: = /, we have 10 possible walks covering the 
(c) For X = c, the weights of walks covering the two-chain subgraph in Fig. [Sff). Since each walks has a weight equal to 
graph are di, dj, dk- Counting the multiplicities of each type — 1^ we have that 
of walk, we have that 

1 V- ^5/ = - E -10 



n 



n 






n ^-^ V 2 / n 

i — 1 ^ ' 2—1 J — 1 



n 

i=l 



where in the last equality we take into account that the number 
where we have used that Y.v,r^v,r^vy, dt = ELi (t)^* ™d of subgraphs of the type depicted in Fig. [^f) and centered at 
'l2v'^vi'^vk '^i ~ 2 Ei=i Ej=i ^ij ('^i — 1) ^j- Thus, node Vi is equal to the number of triangles touching node vi, 

ti, multiplied by {di — 2) (where we have subtracted —2 to 

Tsc = -E^?--E^? + - E a»jrf*rfj--E^' 

2—1 z— 1 l<2,j/<n i — 1 



5 Y^ ,3 5 v^ 72 5 Y^ J J 5 V^ j2 the degree to discount the two edges touching Vi that are part 

n -^ ' n -^ * n .■^^ *"' * ^ n ^ •' of each triangle counted in ti). Hence, we have that 



^ (5,3 - 2^2) + lOQd ^_ _ _10^^^,^10 



T5/ = — y d,u + — }2u 



n -^ — ' n 

4 = 1 j = l 



{d) For X = d, we can determine all possible closed walks _ _, ^^ , pj-, * ; 

of length 5 using the edge graph in Fig. l6td) and derive that 

1 ■r-^ / 2 2 2 (s) For X = /, we have 10 possible walks on the pentagon 



Tsd = - Y. -30 - 10 (d^ + d2 + d; 



Vi^Vj^Vk^Vi 



and the associated weight of each walk is —1. Hence 



+d,dj + djdk + dkdi) 1 ,^ 

on A -in ^5g = - V -10 = -10P/?1, 

^^^ ^° y (3d? + 3d.d,) " ^ 



n 



n n 

n n n 



. ,. ,. ,. where P is the total number of pentagons in Q. 
oOZa oO ^ — A ^ — T ^ — A 2 

2^2^ 2^ iijkd^ Finally, since m^ [Lg) = Tza+T5b + ---+T<ig, we obtain the 

i=i j=i fe=i expression for the fifth Laplacian spectral moment in the state- 

5 JL^ JL^ ment of the theorem after simple algebraic simplifications. ■ 



n n 



y y a„&y(d,dj), 



n 

i=i 3=1 
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(^ ) 



{a) (b) (c) (d) 

Fig. 6. Collection of possible graphs covered by closed walks of length 5. 




"j (d,j^k ^j{^) (djv, -1 v^d/) ri\ -' '7;) 



ie) 
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